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' Abstract. It is shown how data on the cluster correlation function can be used in order to 
reconstruct the density of the pregalactic density field on the cluster mass scale. The method is 
applied to the data on the cluster correlation amplitude - richness dependence. The spectrum 

■ of the recovered density field has the same shape as the density field derived from data on the 
. galaxy correlation function which is measured as function of linear scales. Matching the two 
I amplitudes relates the mass to the comoving scale it contains and thereby leads to a direct 
' determination of Q. The resultant density parameter turns out to be $7=0.25. 
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1. Introduction 



This paper presents another way of determining by comparing density fields determined 
from two independent datasets: APM data on the galaxy two-point correlation function and the 
data on the cluster correlation function - richness dependence. For reasons that will become 
I clear later in the paper, I will reconstruct the quantity which is uniquely related to the correlation 
^ ' function, £,{r), of the density field, or its Fourier transform - the power spectrum P{k): 

>' ^ ^ . 2 

: A2 = {{6M/M f) = Att i{r')r'^dr' = 47r P{k) (^^^|^) k^dk (1) 



I will limit the discussion and the formalism to scales where the present density field is linear 
and therefore can be assumed to reflect the initial conditions, i.e. r > rg = 8/i~^Mpc. 

The outline and the idea of the paper are as follows: in Sec. 2 I discuss reconstruction of (1) 
as function of the linear scale, r, from the APM data. In Sec. 3 I show how to use the data on 
the cluster correlation amplitude - richness dependence in order to reconstruct (1) as function 
of the cluster mass. Comparison of the two density fields gives the amount of mass contained in 
the given linear scale and provides a direct determination of I will show that the two density 
fields, although recovered in different and independent ways, give consistent results and require 
$7 ~ 0.25 with a very small uncertainty. For more details the readers are referred to 

2. Density field from galaxy clustering 

First, I establish the rms fluctuation, A(r), as function of the linear scale r from the APM 
data on the projected angular correlation function w{6). For this I use the APM data on w{9) [ 
0] divided into six narrow magnitude bins Am6~0.5. Galaxies located in each of the bins span 
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a narrow(er) range of so with the data presented in this way one can isolate effects of possible 
galaxy evolution in the entire APM catalog spanning 17<6j<20 or 0.07<2;<0.2 [ ^. 

The projected angular correlation function for each bin is related to the 3-dimensional power 
spectrum of galaxy clustering, P{k), via the Limber equation: 

w{e) =7r r dz{cdt/dz)-^(l){z) r dkP{k; z)kJo{^^^^) (2) 
Jo Jo I + z 

where x{z) is the comoving distance, t is the cosmic time, and the selection function is (j){z) = 
(dN/dz/Ntot)'^. The latter is related to the range of magnitudes for galaxies in each bin, [m^, 
and their luminosity function, via dN/dz = (dV/dz) /^(^^'j ^{L;z)dL with V being the 

comoving volume. Eq.(2) allows to relate the data on w{6) to the underlying power spectrum 
for each of the six narrow magnitude width bins once all of the following are specified: the 
luminosity function in a given band (blue for APM), the i^-correction which accounts for the 
shift along the galactic energy spectrum resulting from cosmic expansion, and the (possible) 
galaxy evolution out to the edge of the APM sample. For the numbers in the remainder of this 
section the luminosity function was adopted from the measurements of [ the i^-correction 
was modeled from standard spectra of galaxy populations [ ^] . 

In order to estimate the importance of the (possible) evolution effects, I proceeded as follows: 
On very small angular scales, the data show that the angular correlation function can be de- 
scribed as a power law, w{6) = Aw6~'^ with 7 = 0.7. On the other hand, as — > 0, the main 
contribution in the Limber equation comes from very small linear scales where the spatial galaxy 
correlation function can be approximated as ^(r)=(r/r^)~^~'>' with = 5.5/i~^Mpc. Relating 
the small-scale ^(r) to w{9) via the Limber equation leads to the following expression for A^- 

where ^{z) accounts for the evolution galaxy clustering, e.g. ^'^(z) oc (1 + z)~^ for clustering 
stable in comoving coordinates. Comparing the value of = w{9)6^''^ computed from eq.(3) 
with the data for each bin allows to constrain the extent to which galaxy evolution affects 
inversion of the APM data in terms of the underlying spectrum. The fits of the power-law 
galaxy correlation function to the small-scale data on w{9) in all six narrow magnitude bins 
show that no galaxy evolution corrections are needed beyond the normal K-corrections and 
evolution of the clustering pattern with time [ |l| . 

Analysis of the data in narrow magnitude bins without galaxy evolution shows that the power 
spectrum obtained by deprojection of the entire APM dataset [ ^ fits the data well in all mag- 
nitude bins. Neglecting the dependence on fibaryon, CDM power spectra can be parameterized 
by only two parameters: the primordial power index n and the excess power parameter VLh [ 
^. Fits of the CDM models to APM data in narrow magnitude bins at all depths show that 
the models require i7/i=0.2 if n=l, or i7/i=0.3 for tilted models with n=0.7. The latter thus 
still requires low in order to fit the APM data in narrow magnitude bins. CDM models with 
larger values of Qh would give smaller w{9) at large 9; whereas smaller values would overshoot 
the data [ [l| . 

Left panel in Fig.l shows the rms fluctuation from the flts to the APM data vs the comoving 



x{z) 



^^{z){l + zfVl + ^zdz 



(3) 
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scale. It is plotted in units of the fluctuation at rs=8h ^Mpc, A(r)/A8. The reason for plotting 
the ratio is that for linear biasing the vertical axis in Fig.l is independent of the bias factor. 

3. Density field from cluster correlation amplitude — richness dependence 

In this section I show how, using the data on the cluster correlation amplitude - richness 
dependence, one can reconstruct the quantity plotted in the left panel of Fig.l as function of the 
cluster mass. Comparing the results with A over a given range of r gives the amount of mass 
in the given comoving scale and leads to a direct determination of O. 

I assume that clusters of galaxies formed by gravitational clustering [ ^ and evaluate their 
correlation function based solely on this assumption [ |9|, |l^. I.e., I assume that clusters of 
galaxies are identified with regions that at some early epoch Zi had initial overdensity such that 
they would turn-around in less than the age of the Universe. It is convenient to choose Zi to be 
sufficiently high when the density field on all relevant scales is in the linear regime. In that case, 
the amplitude, 5ta, of the fluctuation at Zi needed for it to turn around today at z=0 is related 
to A8,i, the amplitude which grows to A8=l at z=0, via 5ta = Qta^%,i- The factor Qta—^-^^ 
for a spherical model and is almost independent of either Zi or cosmological parameters. 

I further assume that the initial density field was Gaussian. In that case the joint probability 
density to find density contrasts 81^2 on scales containing masses Mi^2 respectively is given by: 



p(<Ji; -52) = J J exp{-iq • S) exp(--g ■ C ■ q)d q (4) 

The correlation matrix, C, is related to the spectrum of the primordial density field: its diagonal 
elements are the mean square fluctuation A^ on scale containing mass M and the non-diagonal 
elements are — ^(r) at cluster separations r greater than the comoving scale subtended by the 
cluster masses (~l-3 /i~^Mpc). The probability of two such fluctuations to turn-around by now 
and thereby form clusters of galaxies is PmiM2=Is^^ Js^^Pi^iT ^2)d5id52- The fraction of such 
pairs at the present time would be fMiM2=d'^ PM1M2 / dMi / dM2. The probability for a single 
cluster to form by now is Pm = Js^^p{8)d5] the fraction of such clusters is fM=dPM/dM. 

Now one can construct the correlation function between the present-day clusters of different 
masses. By definition, the 2-point correlation function of an ensemble of objects with number 
density n is given by the probability to find two objects in small volumes dVi,dV2 as dT'12 = 
n^(l + ^)dVidV2. Since the clusters of mass Mi, M2 will make the fraction /M1M2 of such pairs, 
the probability to find them is dVMiNh = /AfiAf2^^i2- On the other hand, the mean number 
density of clusters of mass M would be Jm x n and by definition the probability to find two 
clusters is dVMihh = f&hfAh^'^i^ + (,MiM2)dVidV2. Hence the correlation function of clusters 
of mass M is given by: 

l + ^^HM2 = T^il + ^^) (5) 

J Ml JM2 

Using expansions in terms of Hermite polynomials [ |l^ leads to the following expression 
[0: 

^MM{r) = AMi^mr) (6) 
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where the amphfication factor is: 

oo 
m=0 



with 

Cm, (X) 



^2m 



(^) _|_ -f^m+2(2^) 



(8) 



Here Hm{x) are Hermite polynomials and C=Qta a^m) ' 'which is directly related to the spectrum 
of the primordial density field on scale M. Hence, the data on the cluster correlation amplitude 



- mass dependence [ 13| can be used to invert eqs.(5)-(7) to directly obtain A(M) for 
the primordial density field. The accuracy of these expressions has been confirmed in recent 
numerical experiments [ . 

Since cluster correlation amplitude is known to depend on the cluster richness, M, one has 
to translate the richness into cluster mass. I assume that the two are proportional with the 
coefficient of proportionality normalized to the data on the Coma cluster: 

M(AA) = McomaT^ = 1.45 x lO^^ (^) h-^M^ (9) 

A/Coma \106/ 



Here, I adopted the values for the Coma cluster and richness following [ 14|. This assumes that 
there are no systematic variations in the galaxy luminosity function in clusters of various masses. 

The data on the cluster correlation function can now be used in conjunction with this mech- 
anism for cluster formation in order to set further constraints on the cosmological models via 
eqs.(5)-(7). Indeed, the dependence of the cluster correlation function at a given richness/mass 
on r constrains the power spectrum on the scale of cluster separation, whereas the dependence 
of the cluster correlation amplitude at a given r on the cluster richness/mass constrains the 
shape of the primordial power spectrum on the scale containing that mass. Eqs.(5)-(7) can be 
used in two ways: On the one hand, one can test cosmological paradigms, such as CDM, by 
evaluating the cluster correlation parameters by assuming their power spectrum and comparing 
computed numbers for both (,MM{r) vs r and Am vs M with observational data. Alternatively, 
one can use these expressions in conjunction with data on the cluster correlation amplitude in 
order to directly obtain the spectrum of the primordial density field on scale M independently 
of the assumed cosmological prejudices. The latter, as I show, can also be used to determine 

In terms of the CDM models, I find that no CDM model, whose only free parameters are 
(r2,/i, n), can simultaneously fit three sets of data: the APM data on the galaxy correlation, 
the data on the slope of the cluster correlation function with scale at a given richness, and the 
data on the dependence of the cluster correlation amplitude at a given separation on the cluster 
richness. Namely, the values of (17, h, n) required by fits to the APM data {^Ih ~ 0.2 if n = 1 
and Qh ~ 0.3 if n = 0.7) would be very different from those required by fits to either of the 
other two datasets and vice versa [ |l]. Accounting for a (slight) dependence on ^ibaryon in the 
CDM transfer function [ would not change these conclusions. 

In order to invert eq.(6) in terms of x, and subsequently A(M), I used the data on the cluster 
correlation amplitude - richness dependence from Fig. 2 of [ Il3|. The amplification coefficient 
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for a given richness/mass in eqs.(5)-(7) depends on the value of .^(r) on the scale where it is 
evaluated. The underlying correlation function ^(r) enters eq.(6) via the first and higher order 
terms and can contribute up to ~10-30% to the total amplification. I chose r=25/i~^Mpc on 
which to evaluate the amplification factor from (7). This scale is sufficiently large compared 
to rg to ensure validity of the analysis, but where at the same time ^(r) can be determined 
sufficiently accurately. In the discussion below I adopted ^(25/i~^Mpc)=0.07 in agreement with 
the APM data (cf. [|l|), but the numb ers that follow are not very sensitive to varying the value 
of ^(25/i~^Mpc) within reasonable limits. 




Figure 1. (Left) RMS density fluctuation from APM data is plotted vs the comoving scale. It 
is shown in units of density fluctuation over radius rg. Solid lines correspond to empirical fit 
from [|^; dotted lines to deprojected spectra from [ B]. Three lines of each type correspond to 
1-0" uncertainty in [ ^ and a similar uncertainty in [|3[. Thick dashed line corresponds to CDM 
models that fit APM data in narrow magnitude bins: 0/1=0.2 with n=l or O,h=0.3 with n=0.7. 
(Middle) Plus signs correspond to the primordial density field inverted from data on clus- 
ter correlation amplitude - richness dependence. Solid line shows the fit to the points: 
A(M)=Ag(M/1.7 X lOi^/i-^M0)-O-275. 

(Right) The values of A(r) from the cluster data for three different values of Q: plus signs cor- 
respond to $1=0.1, triangles to $7=1 and diamonds correspond to $7=0.25. Square denotes the 
(by definition) value of unity of A(M)/A8 at Sh^^Mpc. The lines represent density fluctuations 
derived from the APM data redrawn from the left panel in the same notation. 



Once the data on the rhs of eq.(7) for Am vs richness are specified, eqs.(5)-(7) can be solved 
numerically in order to obtain A(M). The middle panel in Fig.l plots the results of this inversion. 
The top horizontal axis plots the values of M at which A(M)/Ag has been evaluated. The 
bottom horizontal axis shows the mass computed according to the normalization to Coma. 
I emphasize again that this method gives directly the pregalactic spectrum as it was at zi 
independently of the later gravitational or other effects. The plot in the middle panel of Fig.l 
shows a clearly defined slope of A(M) oc Af"'^'^''^ corresponding to the spectral index of n ~ 
—1.3. This slope is consistent with the APM implied power spectrum index of w{9) oc 9"^'"^ . 
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4. Determining from the two density fields 

The recovered density field allows one to relate the mass of the cluster to its comoving scale 
thereby directly determining fl. The direct fit to the points in the middle panel of Fig.l gives 

A(M)=A8(M/M8)-° (10) 

with a=0.275 and Ms=1.7 x lO^'^h'^MQ. This fit is plotted with solid line in the middle panel. 
On the other hand, the mass contained in comoving radius rs is M(r8)=6.1 x lO^^Qh^^MQ. 
Equating this with Mg leads to 17=0.28. 

Furthermore, one can determine 0, by comparing A(r) from the galaxy correlation data over 
the entire range of the relevant r with A(M) derived from the cluster correlation amplitude 
- richness dependence. In order to do this I converted the numbers for A(M) to those at a 
given r using that the mass contained in a given comoving radius in the Universe with density 
parameter Q is M(r)=1.2 x 10^^{r/lh~^M-pc)^Qh~^MQ. The right panel in Fig.l shows the 
values of A(r) from the cluster data for three different values of J7: 0.1 (pluses), 0.25 (diamonds) 
and 1 (triangles). The square denotes the (by definition) value of unity of A(M) / Ag at 8/i~^Mpc. 
The lines represent the APM data redrawn from the left panel. One can see that the spectral 
shape of the density field derived from the cluster data is in good agreement with that of the 
APM. The amplitudes of the two fields would match at all r only for ^=0.25. 

5. Summary and conclusions 

In this presentation I discussed reconstructing density field from the cluster correlation am- 
plitude - richness dependence. It was shown that the data can be inverted to obtain the rms 
density fluctuation in the pregalactic density field on scales containing the mass of the clusters. 
The derived density field has the same spectral shape as the density field derived as function 
of the comoving scale from the APM data. Comparing the two amplitudes fixes the amount of 
mass in given comoving scale and allows to determine Q. The value derived from application of 
the method to the data is 0=0.25. This value of O, obtained after normalizing the mass-richness 
relation to Coma, is in good agreement with that implied by the dynamics of the same Coma 
cluster. This further argues that galaxies trace the overall mass distribution in the Universe. 
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